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Abstract 

Formulating the statistical mechanics for a scalar field with non-minimal ^RcfP' 
coupling in a black hole background we propose modification of the original 't Hooft 
"brick wall" prescription. Instead of the Dirichlet condition we suggest some scatter- 
ing ansatz for the field functions at the horizon. This modifies the energy spectrum 
of the system and allows one to obtain the statistical entropy dependent on the 
non-minimal coupling. For ^ < the entropy renormalizes the classical Bekenstein- 
Hawking entropy in the correct way and agrees with the result previously obtained 
within the conical singularity approach. For a positive however, the results differ. 
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1 Introduction 



A common hope during last years is that the study of quantum entropy of black hole can 
shed light on the problem of obtaining a statistical meaning of the Bekenstein-Hawking 
entropy [|1|, 0. A number of different approaches was proposed (see reviews in 0). 
According to 't Hooft |^ (see also ||5|) the entropy arises from a thermal bath of quantum 
fields propagating just outside the horizon. This is a purely statistical calculation treating 
the quantum thermal bath as a system characterized by some energy spectrum. Being in 
equilibrium at a temperature T = (3~^, the states of the system are distributed according 
to Gibbs. An important feature of this calculation is that the density of states becomes 
infinite approaching the horizon. That is why 't Hooft introduced a "brick wall" , a fixed 
boundary staying at a small distance e from the horizon. Assuming that the quantum 
fields do not propagate within this wall, 't Hooft imposed the Dirichlet condition at the 
boundary. A reformulation of this model was suggested in ||^ by using the Pauli-Villars 
regularization scheme. Introducing a number of fictitious fields (regulators) of different 
statistics and masses (dependent on UV regulator /i), it was shown that this procedure not 
only regularizes the UV problem in the effective action but also automatically implements 
a cutoff for the entropy calculation allowing one to remove the "brick wall". Thus, the 
"brick wall" can be conveniently considered as a fictitious boundary the role of which 
is just to make the calculation simpler. 

In an alternate approach the entropy arises from entanglement by means of the density 
matrix obtained by tracing over modes of the quantum field propagating inside the horizon 
i, i, 0, p. 

The powerful method to calculate both the classical and quantum entropies of a hole 



is to apply the Euclidean path integral approach |]I2[. For an arbitrary field system it 
entails closing the Euclidean time coordinate with a period (3 = T^^ where T is the 
temperature of the system. This yields a periodicity condition for the quantum fields in 
the path integral. In the black hole case for arbitrary (3 this procedure leads to an effective 
Euclidean manifold which has a conical singularity at the horizon that vanishes for a fixed 
value (3 = (3h- The entropy is calculated by differentiating the corresponding free energy 
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with respect to P and then setting (3 = (3h- This procedure was consistently carried out 
for a static black hole and resulted in obtaining the general UV-divergent structure of the 



entropy |]T3[-|jT6|. 

For a quantum matter minimally coupled to gravity these three methods to calculate 
the entropy lead to the compatible resuhs 0, 0], Q, H, [0> H and predict the 
similar structure of the UV divergences of the entropy. As was proposed by Susskind and 



Uglum [14 1, these divergences are absorbed in the renormalization of Newton's constant 



and, according to |]T6|, in the renormalization of the higher curvature couplings in the 
effective action. So far, however, there has not been agreement between these approaches 
for the quantum entropy due to a matter non-minimally coupled to gravity. This problem 



was considered in Wh, |18|, [ISl, 20 , 21 



The peculiarity of the non-minimal coupling is easily seen in the conical singularity 
approach [|l^. Consider the Euclidean path integral for a scalar field with the non- 
minimal operator — = — □ + ^R. For an Euclidean manifold with a conical singularity 
the scalar curvature behaves at the singularity as ^-function even if the regular part of 
the curvature vanishes. This 5-like potential modifies the spectrum of the operator — 
and the resultant path integral. The similar phenomenon happens for a manifold with 
boundaries [^. In that case the operator — has a 5-like potential concentrated on the 
boundary. In result, the conical singularity method gives rise to the following entropy due 



to the non-minimal quantum scalar matter |18 



>9, = ^As(l-60ci(/i) 

j^R+^^jj.R,X<-'^R,.apn'tntn)^^^^^ , (1.1) 

where Aj^ is the area of the horizon S and {rij, z = 1,2} is a pair of vectors orthogonal 
to S. The constants Ci(yu) and C2(/i) depend on the UV (energy) cutoff /i and for large 
values of /i they grow to leading order as /z^ and In /i^ respectively. 

The divergences in ( |1 . 1| ) take the correct form so that in combination with the bare 
entropy they are absorbed in the renormalization of Newton's constant and the quadratic- 
curvature couplings in the effective action. In particular, Newton's constant is renormal- 
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ized as follows 123 



7^ = 7^ + (l-60ci(/^) . (1.2) 

One can see that even in flat space there still exists dependence on ^ in the quantum 
entropy (|1 . 1| ) and in the renormalization of Newton's constant (|1.2| ). 

There are simple arguments showing that the quantum entropy in the 't Hooft ap- 
proach has a different dependence on ^. Consider a black hole background satisfying the 
constraint R = 0. Then, solutions of the field equation of motion (□ — S,R)4> = are 
the same as for = 0. Consequently, the energy spectrum of the system and the entropy 
are not effected by the non-minimal coupling. The similar arguments are applicable for 
the entanglement entropy as well to argue that this entropy also does not depend on 
^. This allowed some authors to make a conclusion that "entropy of a quantum matter 
non-minimally coupled to gravity does not have a statistical meaning" . An important cir- 
cumstance missed in this sort of arguments is that a non-minimal matter possesses some 
non-trivial interaction with the horizon (the importance of this interaction was argued 
in |21[]). The "brick wall" prescription is not appropriate to describe this interaction. 



The goal of this paper is to modify the 't Hooft approach in order to get the correct 
^-dependence for the statistical entropy and obtain the correspondence with the conical 
singularity method. 

An important point in our consideration is an idea that in the non-minimal case 
we are obliged to impose certain boundary condition the form of which directly follows 
from the form of the non-minimal coupling. Therefore, our strategy is to replace the 
Dirichlet boundary condition appearing in the original 't Hooft calculation by some ^- 
dependent condition. A motivation for this, actually, comes from the Euclidean version of 
the theory. Indeed, the presence in the operator —D^ of the (5-like potential concentrated 
on the boundary or at the conical singularity can be precisely re-formulated as imposing 
an appropriate condition on the field functions on the boundary or at the singularity 
respectively. In the Euclidean theory these conditions are found to take a simple form 
||24|| . To make our condition in the Minkowskian space-time more clear, note that the 
standard Dirichlet condition for a wave equation means that a wave is reflected by the 
boundary with the change of phase equal to tt. Analogously, what we propose is essentially 
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an ansatz for the field function near the horizon that describes scattering by the hole with 
some non-trivial change of phase. As a result of the scattering, in the spectrum of the 
system appear some low-energy modes the density of which, being proportional to ^, 
grows as g^^^'^{e) approaching the horizon. Namely the contribution of these modes to 
the statistical entropy is essential to obtain the correct dependence on ^ in agreement with 
( |1 . 1|) . In principle, our calculation is quite similar to the standard statistical consideration 
of the non-ideal gas p5[: the corrections to the thermodynamical quantities of the ideal 



gas are expressed via the two-particle scattering phase shift. The only difference from 
our case is that in our model a particle interacts with the horizon rather than with other 
particles. It is important to note that we propose the scattering not just for needs of the 
entropy calculation. We believe that it is an actual interaction between the non-minimal 
matter and the horizon that is dictated by the form of the non-minimal coupling. 

2 WKB field function and boundary condition 

A straightforward generalization of the approach [Q, |P] is to consider the more general 
case of static black hole with the spheri-symmetric metric written in the form: 

ds^ = -gir)de + -^dr^ + r^idO^ + sin^ Od^^). (2.1) 

Not assuming this metric to satisfy any concrete field equations, we suppose that g{r) 
is a non-negative function having a simple zero at r = r+ where it behaves as g{r) = 
;^(r — r_|_) + 0((r — r_|_)^). This corresponds to the position of the outer event horizon 
at r = r+. The Hawking temperature calculated for the metric (|2.1|) is Th = Pjj^- Note 
that the scalar curvature R for the metric ( p.l| ) is function of only radial coordinate r. 

In this black hole background we consider a non-minimal scalar field which satisfies 
the equation: 

(a - - ^R)(f) = . (2.2) 

Expanding the wave function in spherical coordinates (f) = e^^^Yim{0 , ^p) f (r) , we obtain 
the equation on the radial function 



gV/(r) + r Vr(r Vr/(r)) - r^g{m^ + ^R+ ^^^^ )/(r) = . (2.3) 
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The WKB approximation provides us with the following solution of the Eg. ( |2.3| ): 



/(r) = p(r) e^'f T'^i-^ , k{r) = ^E^-{m^ + iR + ^-^^±^)g{r) , (2.4) 

which is valid in the region where k'^ir) > 0. It is clear that for large mass m this is space 
closely located to the horizon. If R takes a non-zero value at the horizon then this region 
can be approximated by some sort of the constant curvature space. However, the exact 



results for such a space show [^] that the mass parameter in the solution enters only in 
the combination (m^ — \R)- The WKB approximation, being well-defined for large values 
of /, does not give this result. Correcting this, we arrive at the following wave function 



/(r) = p(r) (e^ff'^M + Ae-^rf^M 



k{r) = p^-{M\r) + ^-^^^)g{r) , (2.5) 

where we denote M^(r) = m? — — ^)R. Since the quantity M^(r) is slowly varying 
near the horizon it can be considered constant there. The constant A in (|2.5|) is to be 
determined from the boundary condition. The amplitude p(r) is a slowly varying function 
of r. In what follows, we neglect its derivatives with respect to r and omit writing it in 
the formulas. 

Consider a boundary staying at a small distance e from the horizon S. In the limit 
e — * it approaches the horizon ^ S. The parameter e is assumed to be smalleiQ 
than any UV cutoff fi~^ appearing under regularization of the quantum field theory 
and such quantities as /^^5'(e) are considered to be negligible. Therefore, the boundary 
plays an intermediate role just to simplify the consideration. The condition, in principle, 
can be imposed directly on the horizon. 

In order to arrive at an idea of the boundary condition to be imposed on let us 
start with the following "simple-minded" condition: 

(n^a^0-e«:0)|s. = , (2.6) 

where n'^ is vector normal to and k is the extrinsic curvature of E^, k = V^n'^. Namely 
the condition of this type appears in the Euclidean version of theory p4| . For the metric 



'''Explicitly, this means that << 1. 



( p.l| ) and sufficiently small e the condition ( p.6| ) reads 

((7(r)9J(r) - Cf{r)) |.=.^+. = , (2.7) 

where ^* = 27i(3jj^^. For the function ( |2.5D it gives the expression for the constant A: 
A = e"'*^'^-'. The phase r](k) is defined as follows: 

where /i;(e) = k{r = r+ + e). 

The condition we are looking for is a ^-dependent modification of the Eq. (|2.7|) . Curi- 
ously enough, an appropriate condition takes the form 

(/7(r)a,-rr/(r)U++e = 0,z/ = 2^ , (2.9) 

where {g{r)dr — is a pseudo-differential operator. It acts as follows 

where we neglected the derivatives of k{r). The condition (|2.9| ) leads to the following 
constant 

A = e"'"'''''^ e*^(''+^) , (2.10) 

where ri{k) = arctan (^J^^^r?) is defined as in ( p.8| ). Remarkably, for z/ = (^ = 0) the 
Eq. (|2.9|) coincides with the Dirichlet condition. 

With the condition (|2.9|) imposed, the WKB field function takes the form 

/(r) = (^e^^++' + e*^''('=)e*^(^+^)e~*^'-++= . (2.11) 

This is that ansatz which we propose for the field function near the horizon. The function 
( p.ll|) is valid in the region r+ + e < r < te-, where te is determined by the equation 
k'^i'^E) = 0. The precise meaning of the condition appears when we consider the field 
function at the horizon. Introduce new coordinate z = f g~^dr. Then, at the horizon {z 
goes to — oo) the function ( p.ll| ) becomes 



■lE - C 



where zo is a constant, and describes scattering by the hole with change of phase i/ri{E) + 
7r(/y + 1). 



7 



3 Density of states and entropy calculation 



In order to discretize the energy spectrum and simplify counting of states we impose also 
the Dirichlet condition = for r = te- This gives the quantization condition 

f-B dr 

2 —-kir) = vr]ik) + 2-nn + -Kv , (3.1) 
Jr++e g[r) 

where n is an integer number. Inverting this, we get the equation 

n = n{E,l)=no{E,l) + ni{E,l) , 
no{E,l) = - / —k{r) , 

= (arctan ( ^(^^^.2 ) + ^) , (3-2) 

solution of which gives the energy levels {En,i} at a fixed /. 

Near the horizon the metric function can be approximated a.s g{r) ~ ^^(r — r+). In 
this approximation the integral in ( p.2|) is calculated exactly. We find in result that 

n„(E,l) ^'^U arctanh^l - - {E^^,„ 

„.(E,0^-^(a.cta„(e^^5^J|.)+.) , (3,3) 

where E'La , = (M^ + ''^^t^^ )g(e). We see that E = Emini is the minimal possible energy 
at a fixed /. For E ^ Emin,i the behavior of the function n{E, I) (|3.2|) is mainly determined 
by the component ni{E, I) while for large enough E the component no{E, I) becomes the 
leading one. For 1/ = 2^ the n{E, I) is monotonically increasing (both for ^ > and ^ < 0) 
function]] taking value n{E, I) = ^ for E = Emin,i- The total number of modes with energy 
less than E can be determined by taking the quantity {n{E, I) — n{Emin,i)) and summing 
over the degeneracy of the angular modes: 

n{E) = J dl{2l + l){n{E, I) - n{E^,^j)) = no{E) + n^{E) , (3.4) 

^ In general the behavior of the function ni{E,l) depends on relative sign of v and ^. Namely, if 
sign{u) = sign(^) the function ni(E,l) is increasing for E > Emin,i while in opposite case, sign{v) = 
sign{^), it is decreasing. This affects the behavior of the whole function n{E,l) which is monotonically 
increasing in one case {sign{u) = sign(S.)) and has one minimum in other case {sign{y) = —sign{^)). 
Here we assume that = 2^ and discuss the possible choice u = — 2|^| below in the Discussion section. 
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where the sum over / has been approximated by an integral, and this integration runs 
over non-negative values of / for which the square roots k{r) and k{e) in the integrand are 
real. After the integration we have 



no{E) = - / dr^(E' - M'g{r)f' (3.5) 



and 



V r 



2 



ni{E) = / dkk aictan 



TT g{t) Jo \k'^ — ^ 



.2 



'-^{^ arctan + i*k{t) - ^ arctan , (3.6) 



where fc(e) = _ M'^g{e). 

One can think about this system of particles as consisting of two components: the 
ordinary particles with the number rio(-E') as in the original 't Hooft model and the 
scattering particles with the number n\[E). Remarkably, n\{E) is proportional to the 
horizon area ^4+ = 47rr^. 

For E ~ Emin = Mg^^'^(e) the functions uq^E) and ni(E) behaves as follows 

no{E) ^^rl^iE^ - ElJ'/^ 
15 vr h. 



mm 

n,{E) ^ szgn{u0^rl^{E' - E^J^' . (3.7) 

min 



If sign{v^) = +1 the total number of modes n{E) is the monotonically increasing function 
of E and E = Emm is actual minimum of energy. 

To determine the thermodynamics of this system, we consider the free energy of a 
thermal ensemble of scalar particles with an inverse temperature (3 

PF = / dE^ ln(l - e-^^) = PFo + PF, , (3.8) 

^ ^min (111/ 

where we separate the contributions due to the ordinary modes and the scattering modes. 

Applying Pauli-Villars regularization scheme for the present four- dimensional scalar 
field theory, one introduces five regulator fields {0,, z = 1, 5} of different statistics and 
masses {mj, i = 1,...,5} dependent on the UV cut-off fi 0. Together with the original 
scalar 0o = ("^0 = fn) these fields satisfy two constraints: I]i=o = and J2t=o ^i'l^f = 



0, where Aj = +1 for the commuting fields, and A, = —1 for the anticommuting fields. 
Additionally, we assume that all the fields have the same non-minimal coupling = ^, i = 
0, ...,5. Not deriving the exact expressions for rrii, we just quote here the asymptotes 

Ajm^ In = fi^bi + In — - + m'^b2 , 

^Ailnm^ = \n— , (3.9) 
j=o 

where bi and 62 are some constants, valid in the limit — > 00. With contribution of each 
field added the free energy ( p.8|) becomes 

_ 5 

PF = Y,A,f3F' = pFo + (3F, . (3.10) 

i=0 

For the scattering modes of a single field we have 



PF, = r dE^ ln(l - e- 

^ ^min (111/ 



= A f dEEarctan ( ^^^^^^^ 1 1^(1 " • O-H) 

A remarkable property of the expression ( |3.11D is that for v = 2^ it depends only on the 
absolute value of ^. Integrating over E in ( |3.11|) and focusing only on the divergent for 
small e terms, we find 



r 



2 



/3F, ^ -\^\-t-(El^JniPE^,^) - i-ElJ - 1^ , (3.12) 
^(e) 2 p^g{e) 

where x is an independent on M constant. Summing contributions of each field (pi, we 
obtain 

PF,^-^-frlJ2^^Mf\nM^ , (3.13) 

where = - (| - 

Our calculation of a part of the free energy due to the ordinary modes essentially re- 
peats for more general static metric ( p.l|) the calculation presented in 0. In the expression 
for the free energy of a single field we can integrate by parts and get 

pF, = -(3r ^W^dE . (3.14) 
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Focusing on the divergences at the horizon in the expression for the number of ordinary 
modes no{E) (P75|) we find 



(3.15) 

where C = {R^yn^n''^ — 2R^yapn^nfn^jn^j)\Y. with notions as in ( p..l| ). 

Substitution of this in ( |3.14|) , integration over E and summation for all fields 0j give 
us the total free energy due to the ordinary modes: 

5 1 1 5 



^ -if .igA,M^nM^ - A.,.>A^ . (3,16) 

It is a manifestation of the mechanism discovered in [^] that the expressions (|3.13| ) 
and ( p.l6|) are regular in the limit e ^ 0. There is a precise cancelation of the divergences 
between the original scalar and regulator fields. The resultant expressions, however, 
become dependent on the UV regulator /i. 

Altogether, ( p.l3|) and (|3.16|) give the total free energy ( p.lO| ) of the system. Calcula- 
tion of the entropy S = at the Hawking temperature = j3jj^ gives 

S = \{\ + \i\)rlj:^Mf\nMf + ^\rlcj:^dnMf . (3.17) 

Using the definitions of C and M'^ = — (| — C.)R, and assuming that the value of 
the scalar curvature R at the horizon is much smaller than each mf, we arrive at the 
expression 

s = Ia. ^(1 + m)±^.m^ m™? + {-^(i + ml -of^R 

+^3^/^(fi„X<-2B„„„a<<n>f)}i:AihimJ . (3,18) 

This is our main result. In minimal case (^ = 0) and for the Reissner-Nordstrem back- 
ground Eq.( p.l8| ) coincides with the result of p. 



4 Discussion of the result 



It should be noted that the expression ( p.l8| ) is not exactly that we were going to obtain 
for the entropy anticipating the complete agreement with the result (|1.1|) found within the 
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conical singularity method. The dramatical difference of our result is that the statistical 
entropy ( |3.18| ) is not analytic with respect to the non-minimal coupling ^. For negative ^ 
it exactly reproduces the conical expression (|T]l|) while for positive ^ the both quantities 



differ. This is well illustrated in the leading order. The conical entropy (|1 . 1|) then takes 
the form Sq = ^A-^^l — 6^)ci(/i) and becomes negative for ^ > |- This puzzling behavior 
has been discussed in and more recently in [|^. In particular, it was noted 

that a statistical entropy defined as S* = —Trp In p for a density matrix p can not behave 
in this way being automatically positive. Therefore, may be it is not so surprising that 
our semi-classical computation of the statistical entropy gives rise to the expression 5" = 
+ 6|(^|)ci(/i) which is always positive. 

However, can we modify in some way our calculation and reach the complete agreement 
with dLlI)? At first sight, it seems possible if we assume that u = —2\^\ instead of u = 2C, 
. Then the free energy f3Fi looks becoming dependent on sign of ^ and the entropy 

does too. However, in this case the function n{E,l) ( p. 21 ), ( p. 3D is not monotonically 
increasing for ^ > 0. Instead, it develops a minimum in a point Emin,i different than 
Emin,i- One observes the same behavior for the function n{E) ( p.4|) having sense of the 
total number of modes with energy less than E. As is seen from ( p.7| ), for sign{i/^) = — 1 
it decreases at E = Emin, takes minimal value at some point Emin and then monotonically 
increases. This behavior means that for u = —2\^\, ^ > we have to re-count the number 
of modes which is no more given by the expression ( |3.4| ). In result, we should get new 
function n{E) which is monotonic and leads for ^ > to a positive entropy as well. 

The computation we present in this work sheds some light on the origin of the ^- 
dependent part of the entropy of the non-minimal scalar field. It is important to note that 
the main contribution to (3Fi comes from modes with energy close to Emin = Mg^^'^(e). 
For this energy we may substitute ln(l — e'^^) ~ \n{j3E) in (|3.11|) and this is the divergence 
(if E —>■ 0) that appears in (|3.12|) . It is interesting that one may give an alternative simple 
calculation of /3Fi assuming existence in the spectrum ( |3.1| ) exactly^ n = |^| modes with 
the minimal energy Ei = Emin,i = (M'^ + ^^^-2^)*-^^^''5'(e). Indeed, summing then over / we 



^ Since n is an integer it is better to suppose that n = [£] where [£] is integer part of ^. For a big ^ 
this is, however, not important. 
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find 



= / dl{2l + l)|e| ln(l - e'^^O = H dEE\n{l - e^^ 



-'min — -Cj; — 



EUn(/3E™,„) (4.1) 



that is exactly the expression ( ^.121) . So, namely, these low-energy modes are responsible 



for the (^-dependent part of the entropy. The role of the boundary condition we propose in 
the Section 2 is just to provide us with necessary number of these modes. Unfortunately, 
it is not clear that we have a well-posed problem with those boundary conditions and, 
for example, the Hamiltonian is obviously self-adjoint (though we deal with apparently 
positive energy spectrum). Therefore, it may happen that there exists a more strict way of 
getting the low-energy behavior we mentioned above within a well-formulated self-adjoint 
extension of the scalar field Hamiltonian. 
The following conclusions are in order. 

1. The quintessence in understanding the entropy of the non- minimal scalar field 
is existence of the low-energy modes number of which is governed by the non-minimal 
coupling ^. In our approach these modes appear due to the non-trivial scattering con- 
dition which we impose on the scalar field at the horizon. In the limit e ^ ("brick 
wall" removed) they presumably become zero modes of the Hamiltonian which are in- 
finitely degenerate due to the angle dependence. However, their status in the well-defined 
Hamiltonian picture is not clear. 

2. For negative ^ the entropy ( |1.1|) may have a statistical explanation within the 
procedure described in this paper. However, we still lack this explanation for (|1 . 1| ) for 
positive ^. Moreover, the validity of the expression ( |1 . 1| ) for ^ > is under question. More 
careful analysis |^ shows that the conical method is not unambiguous and one may be 



needed to re-consider the way one obtains the expression (|1 . 1| ) within this method. 

3. For ^ < the divergences of the statistical entropy may be renormalized by the 
renormalization of Newton's constant according to (|1.2| ). However, for ^ > even after 
such a renormalization we still have a divergence in the entropy behaving as Sdiv = 
3|(^|v4sCi(/i). It is not clear in a renormalization of which quantity it may be absorbed. 

4. Our result may be considered as a confirmation (for the non-minimally coupled 
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scalar field) of the point that any physical quantum field (bosonic or fermionic) must have a 
positive entropy. This, however, may have dramatical consequences for the understanding 
the origin of the black hole entropy in theories of the induced gravity. In this kind of 



theories ^Tj the Einstein gravity arises in the low-energy regime by averaging over the 
constituent matter fields interacting with a background (classical) metric. The set of 
these fields is specially arranged to make the induced Newton's constant Gind UV finite. 
It is hoped p8|, p2| that the finite black hole entropy Sbh = jtt — can be induced 
in a similar way as a statistical entropy of the constituents. This hope is essentially 
based on the possibility to make the black hole entropy finite by making finite Newton's 
constant. However, this assumes that among the constituents there should present fields 
carrying negative divergent entropy which compensates the positive divergent entropy of 



other constituents. In a concrete realization of the induced gravity the role of those 
particles is played by the scalar fields non-minimally coupled to gravity with positive ^. 
However, if each constituent has a positive divergent entropy the total induced entropy 
is also positive and divergent even if the induced Newton's constant is finite. Thus, we 
have an obvious problem with inducing the Bekenstein-Hawking entropy. Possibly, the 
approached developed in can be useful in resolving this problem. 

One certainly needs understand better all these issues and we hope to return to them 
in further publications. 
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